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ABSTRACT 

We present spectral analysis of NuSTAR and Swift observations of Cep X-4 during its outburst in 2014. We 
observed the source once during the peak of the outburst and once during the decay, hnding good agreement in 
the spectral shape bet'ween the observations. We describe the continuum using a po'werla'w 'with a Fermi-Dirac 
cutoff at high energies. Cep X-4 has a very strong cyclotron resonant scattering feature (CRSF) around 30keV. 

A simple absorption-like line 'with a Gaussian optical depth or a pseudo-Lorentzian profile both fail to describe 
the shape of the CRSF accurately, leaving significant deviations at the red side of the line. We characterize 
this asymmetry 'with a second absorption feature around 19keV. The line energy of the CRSF, which is not 
influenced by the addition of this feature, shows a small but signihcant positive luminosity dependence. With 
luminosities between (l-6)xl0^®ergs ', Cep X-4 is below the theoretical limit where such a correlation is 
expected. This behavior is similar to Vela X-1 and we discuss parallels between the two systems. 

Subject headings: accretion, accretion disks — radiation; dynamics — stars; neutron — X-rays; binaries — 

X-rays; individual (Cep X-4) 


1. INTRODUCTION 

Neutron star high-mass X-ray binaries (HMXBs), i.e., neu¬ 
tron stars accreting from an early-type stellar companion, typ¬ 
ically show very high levels of variability in their X-ray emis¬ 
sion. As the unabsorbed flux is directly related to the mass ac¬ 
cretion rate, we can understand these changes as originating 
from variations of the latter. The physical conditions inside 
the accretion column, where most of the X-rays are produced, 
change with accretion rate. By sampling different luminosity 
levels we can obtain constraints on the geometry and physical 
processes in the accretion column. In particular, Be-star sys¬ 
tems, wh ere the stellar companion has a large circumstellar 
disk (e.g. [Okazaki & Negueruela|200T] l, are ideally suited for 
studying the luminosity dependence of the X-ray spectrum, as 
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they show weeks- to months-long outbursts which can cover 
more than two orders of magnitude in luminosity. 

Many accreting neutron stars show prominent cyclotron 
resonant scattering features (CRSFs). CRSFs are produced 
by resonant scattering of photons off electrons moving per¬ 
pendicular to the magnetic held. The electrons are quantized 
on Landau levels which energies depend directly on the lo¬ 
cal magnetic held, and therefore make the observed energy of 


the CRSF a direct tracer of the B-fleld (see, e.g., Schonherr 
jet al. 120071 ^nd references therein). With a variable accretion 
rate the cyclotron line production region can move along the 
accretion column, sampling different magnetic fields. A pre¬ 
cise measurement of the CRSF therefore probes the accretion 
geometry of a source. 

According to theoretical calculations, a strongly asymmet- 


emission wings, (e.g. |Yahel| 19791 

Isenberg, Lamb & Wang 

1998l|Araya & Harding|1996 

1999HAraya-G6chez & Hard- 

ing|2000l Schonherr et al.|2007[ Schwarm|2010 among oth- 


geometry and the underlying X-ray continuum. The asym¬ 
metry originates from photon-spawning, where electrons ex¬ 
cited to higher harmonic Landau levels cascading back to the 
ground stat e and emitting photons close to the fundamental 
line energy (Schonherr et al.|2007|l. 


Observationally, however, only very few of the approxi¬ 
mately 25 known CRSF sources show any evidence for a 
deviation fr om very simply-shape d fundamental lines. In 
GX 301-2, |Kreykenbohm| ( |2004| l found marginal evidence 
for an asymrnetric profile at certain pulse phases using RXTE 
data, but the st atistics did not allow for a deta iled description 
of the feature. |Pottschmidt et al. ( 2005| l and Nakajima, Mi-j 


jhara & Makishima| ( |2010| l showed that the fundamental line 
in V0332-I-53 is better described by two absorption-li ke lines 
at almost the same energy, but with different widths. [Iwakirij 
jet al.j ( |2012| ) claimed the detection of a CRSF in emission 
rather than absorption at certain phases of 4U 1626-67, us¬ 
ing Suzaku data. This feature might lead to a complex line 












































2 


Fiirst et al. 


profile in the phase-averaged data, but their data did not pro¬ 
vide a high enough quality to measure it. 

Among the six CRSF sources studied w ith the Nuclear 
Spect roscopic Telescope Array (NuSTAR, [Harrison et al.| 
2()13|l so far, none show significant deviations from a Gaus- 
sian or Lorentzian optical depth profile (jFiirst et al. 2013[ 


Fiirst et al. 2014a IFiirst et al.|2014bt|Bellm et al.|20l4 |Ten-| 

dulkar et a .[2014 [Bhalerao et al.|2015| l. A detailed study of 

Her X-1 revealed good agreement in the line profile between 
NuSTAR and Suzaku, and put the most stringen t limits on 
possible emission wings to date ([Fiirst et al.|2013|) 

The Be HXMB Cep X-4 has a CRSF around 30 keV, which 
is ideally suited to be studied in detail with NuSTAR, thanks 
to the instruments unprecedented energy resolution as well 
as increased sensitivity above lOkeV compared to p revious 
missions. C ep X-4 was discovered by OSO 7 in 1972 ([Ulmer 


et al. 
et al. 


\912) and again detected by Ginga in 1988 (Makino 


^88[l! During the 1988 outburst, regular pulsations with 


a pulse period around 66 s were discovered, and evidence for a 


CRSF around 30keV was found (|Koyama et aL|1991 

[Mihara 

et al. 1991||. The optical counterpart was identified by 

Bonnet- 

Bidaud & Mouchet|([1998[l, who measured a distance of 3.8 + 


0.6 kpc. 

The most detailed spec tral description to date is presented 
by [McBride et al.[([2007|l, who used Rossi X-ray Timing Ex- 
plorer (RXTE) data taken during an outburst in 2002. They 
confirm the CRSF around 30.7 keV and describe the contin¬ 
uum with an absorbed power-law with a Fermi-Dirac cutoff. 
By monitoring the source over the outburst they find a hard¬ 
ening of the broadband spectrum with luminosity, but the data 
quality does not allow the investi gation of the dependen ce of 
the CRSF energy on luminosity. [McBride et al.| ([2007[l also 
show that the pulse profile changes signihc anTIy as a func¬ 
tion of luminosity, confirming the results by [Mukeijee et al.[ 
POOO I who use RXTE and Indian X-ray Astronomy Experi¬ 
ment (IXAE) data. 

The rest of the letter is organized as follows; in Section]^ 
we detail the data reduction, and in Section]^ we present the 
spectral analysis. Section l^discusses and summarizes our re¬ 
sults and concludes this letter. 

2. OBSERVATIONS AND DATA REDUCTION 

In June 2014 MAXI de tected renewed activity from Cep X- 
4 ( Nakajima et al.|20T4 i, and Swift performed pointed XRT 
observations after an automatic BAT trigger ([Evans et al. 


2014 1 . Ozbey-Arabaci et al. (20141 performed optical ob- 
servations and found evidence for a strong Be-disk, expected 
to occur with the onset of the X-ray activity. We trig¬ 
gered NuSTAR observations and observed Cep X-4 twice, 
on 2014 June 18th-19th (MJD 56826.92-56827.84, ObsID 
80002016002, observation 1) close to the maximum of the 
outburst and on 2014 July lst-2nd (MJD 56839.43-56840.31, 
ObsID 80002016004, observation 2) during the decline. Both 
NuSTAR observations were supported by Swift/XRT snap¬ 
shots (ObsIDs 00080436003 and 00080436004, respectively). 
Figure [T] shows the light curves of Swfft/BAT and MAXI, and 
the average count-rate of the NuSTAR observations. 

2.1. NuSTAR 

We extracted NuSTAR data using the standard 
nupipeline software vl.4.1 as distributed with HEASOFT 
6.16. We used CALDB v20150316, taking a time-dependent 
gain change into account. Standard screening of the data 
resulted in good exposure times of 40.5 ks for observation 1 







Fig. 1. — Lightcurve of the 2014 outburst of Cep X-4 as observed with 
Swift/BAT (blue diamonds), MAXI (green circles) and NuSTAR (red dia¬ 
monds). All rates were rescaled to mCrab in the respective energy band. 
The right hand jc-axis gives the measured count-rates for the NuSTAR data. 
MAXI did not take data during the maximum of the outburst. 

and 41.2 ks for observation 2. We extracted source spectra 
separately for FPMA and FPMB from a circular region with 
a radius of 120" centered at the J2000 coordinates. The 
background spectra were extracted from the opposite corner 
of the NuSTAR field of view from a circular region with 90" 
radius. 

Spectra were modele d with the Interactive Spe ctral Inter¬ 
pretation System (ISIS, Houck & Denicola 1 2000 1 vl.6.2-30 
and errors are given at the 90% conhdence level unless oth¬ 
erwise noted. The data were rebinned within ISIS to a S/N 
of eight between 3.2^5 keV and a S/N of three above that, 
binning at least two channels together. While the NuSTAR 
calibration is nominally good down to 3 keV we choose to ig¬ 
nore the first few bins in order to allow Swfft/XRT to drive the 
model at the lowest energies. We include data up to 60keV, 
where they become background domianted. 

2.2. Swift/XRT 


Dat a from the Swift X-ray Telescope (XRT, Burrows et al. 
[2005 [ l were e xtracted following th e steps outlined in the Swift 
user’s guide ([Capalbi et al.|2005jl using HEASOFT 6.16. Ob¬ 
servation 1 (MJD 56826.94-56826.956) was performed in 
windowed timing (WT) mode, and we extracted the source 
spectrum from a rectangular region with 45" length perpen¬ 
dicular to the read-out direction. The background was ex¬ 
tracted from similar regions on both sides of the source loca¬ 
tion. After standard screening we obtained a good exposure 
time of l.Oks. 

Observation 2 (MJD 56839.462-56839.531) was per¬ 
formed in photon counting (PC) mode and was heavily piled- 
up. Following the procedure described in the XRT data anal¬ 
ysis guid^^we determined that an annulus extraction region 
with inner radius 12" and outer radius 60" removes most pile- 
up while still providing a good S/N spectrum. The back¬ 
ground was estimated from a large region south-east of the 
source. The observation resulted in 1.6 ks of good exposure 
time. Both Swfft/XRT spectra were rebinned to a S/N of six 
throughout the used energy range of 0.8-10 keV. 

3. SPECTRAL ANALYSIS 


As [McBride et al.[([2007|l demonstrate, the hard X-ray con- 
tinuum of Cep X-4 is well described by an absorbed powerlaw 

http://WWW.swift.ac.uk/analysis/xrt/pileup.php 
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Distorted cyclotron line 


TABLE 1 

Parameters of the best-ht Fermi-Dirac cutoff model for 

BOTH OBSERVATIONS. 


Parameter 

Obs. I 

Obs. 2 

Ah [10^2 cm-2] 

1 05+®" 

0061+A?-08 

U.UDI qqjq 

0.83/®;®/ 

1.41 ±0.25 

A 

^cont 

F 

^•^^^-0.005 

0 96+0-^^ 

Ecux (keV) 

24 ±4 

25 ±4 

fifoid (keV) 

-0.6 

C "7+0.6 
-0.8 

£crsf (keV) 


90 4,9+*^-27 

0 24 

o-cRSF (keV) 

5.8 ± 0.4 

4 9 + 04 

(keV) 

20/5 

16 

-3.0 

£ab. (keV) 

19.0/®-5 

18.5 ±0.7 

(Tabs (keV) 

2.5 ± 0.4 

2.1 ±0.5 

rfabs'’ (keV) 

u.ou_o 17 

0.37/®;" 

A(FeKQ-)“ 

(i.39/«|«)x10-3 

(2.8/“) X 10-^ 

o-(Fe Ka) (keV) 

0.42 ± 0.05 

-0.10 

£(Fe KQ')(keV) 

6.474/«;030 

6 39+0.06 
^•-^^-0.07 

/1bb‘= 

(2.22/®") X 10-3 

(7.3/[-3) X lO-'* 

kTm (keV) 

0 899+0.030 

0.96 ± 0.06 

Cfpmb 

1.0319 ±0.0019 

1.023 ±0.004 

Cxrt 

0.962 ± 0.019 

0.91 ± 0.05 

y/d.o.f. 

1215.83/1082 

776.35/670 

D. 

1.124 

1.159 


^ in photons keV ^ s ^ cm ^ at 1 keV 


line depth, optical depth = dl{(T^flk) 
^ in 10^^ erg s~^ for a source at lOkpc 


with a Fermi-Dirac cutoff (|Tanaka|1986]l of the form 


F(E) oc X 1 


(E-E,, 


( 1 ) 


To that continuum, [McBride et al.| (|2007t add a CRSF 
modeled by a multiplicative absorption Im^] with a Gaus¬ 
sian optical depth profile (gabs in XSPEC) and a narrow 
additive fluorescent Fe Ka line. We apply this same model 
to both observations separately, modeling Swfft/XRT, NuS- 
TAR/FPMA, and FPMB simultaneously and allowing for a 
cross-calibration constant relative to FPMA for FPMB and 
XRT (Cfpmb and Cxr t)- We use the phabs absorption m odel 
with abundances by [ Wilms, Allen & McCray 
cross-sections by Verner et al.|p996|k 


(20001 and 


The model results in an unacceptable fit for both observa¬ 
tions (;^f^/dof = 2275/1087 = 2.09 for observation 1 and 
;^f^/dof = 1128/675 = 1.67 for observation 2) with strong 
residuals below 10 keV. We therefore add a blackbody compo¬ 
nent with 1 keV which improves the fit significantly. 

We obtain ;i^^/dof = 1324/1085 = 1.22 for observation 1 and 
— 802/673 = 1.19 for observation 2. The residuals of 
observation 1 for this model are shown in Figure!^. 

As can be seen in Figure 1^, the residuals stillshow some 
structure between 10-20keV which is not modeled by the 
CRSF. We therefore add another multiplicative absorption 
line with a Gaussian optical depth profile. This addition pro¬ 
vides a significant improvement and results in a good fit, with 
X^/dof - 1215/1082 = 1.12 for observation 1 dxvdx^ldoi - 
116I61Q - 1.16 for observation 2. The energy of the CRSF 
does not change significantly when adding the second absorp¬ 
tion model, which is found to be around Eabs ~ 19 keV in both 
observations. The residuals to this best-fit model are shown in 


We note that while a CRSF is produced by resonant scattering, not ah- 
sorption, a possible parametrization is the same as for the latter. 


profile in Cep X-4 

Energy (keV) 



Energy (keV) 

Fig. 2.— (a) Count spectra and best-fit models of both observations of 
Cep X-4. SwiftpiKT data are shown in magenta, NuSTAR/FFMA data in 
red, and NuSTAR/FFMB data in blue. The best-fit model is shown in gray, 
the model evaluated without the absorption lines in green, and the model 
with only the primary CRSF in orange. The inset shows a zoom on the red 
edge of the CRSF in observation 1. (b) residuals of observation 1 to a FD- 
cut model with only one absorption line, (c) residuals of obseiwation 1 to the 
best-fit model, (d) residuals of observation 2 to the best-fit model. Data were 
rebinned for plotting puiposes. 


Figure]^ and d for observation 1 and 2 respectively, and the 
best-fit parameters are given Table [T] 

Using two lines with Gaussian optical depth profiles de¬ 
scribes the shape of the CRSF in both observations very well. 
Using a pseudo-Lorentzian profile instead (modeled by the 
cyclabs model) results in a very similar fit and also requires 
a second absorption feature. The line energies are about 2 keV 
lower, consistent with the expected d ifference between the 
models (see, e.g. [Staubert et al.||20T^ . Our model is there¬ 
fore insensitive to the slight differences in shape betwe en the 
two profiles, similar to results obtained for Her X-1 (jFiirst 
e rar]|20T3ll and V 0332-F35 (jPottschmidt et al.||2005i [feE^ 


jima, Mihara & Makishima|2010|l. 

The CRSF energy shovvs small but statistically signifi¬ 
cant variations with luminosity, decreasing from £’ob.s i = 
30.39 /q[ 4 keV in the first observation to £’obs 2 = 
29.42/Q24keV in the second observation. This behavior is 
qualitatively the same with or without the second absorption 
feature. In order to rule out the possibility that small changes 
in the continuum parameters influence the measured energy of 
the CRSF, we performed a simultaneous fit of both datasets. 
This was possible as the continuum parameters (besides the 
normalization) do not change significantly between the two 
observations. In the simultaneous ht, the photon index F, the 
folding energy Ema, the cutoff energy Ecut, the absorption col¬ 
umn Ah, and the Fe Ka line parameters are tied between both 
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datasets. We obtain very similar results with respect to the two 
absorption features and the black body component, in partic¬ 
ular the luminosity dependence of the CRSF energy is seen 
with the same significance and both observations require a 
second absorption feature. 

In order to further investigate the choice of the continuum 
on the shape and ener gy of the CRSF , we also model the data 
with an NPEX model ( |Mihara| 19951 1, in which we fix the sec¬ 
ondary power-law index to F 2 = -2. We obtain a very sim¬ 
ilar statistical quality of fit and consistent parameters. The 
energy of the line is fitted to Fobs 1 = 30.59^g [jkeV and 
Eobs 2 - 29.43;^q 22 keV. The secondary absorption feature is 
also clearly visible in the residuals and its addition leads to a 
similar improvement in terms of;^^^. 



4. DISCUSSION 

We have presented simultaneous Swift and NuSTAR ob¬ 
servations at two different luminosities taken during an out¬ 
burst of Cep X-4 in June and July 2014. We describe the 
broad-band spectrum with a power-law attenuated by absorp¬ 
tion and a Fermi-Dirac cutoff. We find that the continuum 
does not change significantly between the observations but 
that the CRSF shows complex behavior. 

4.1. Shape of the CRSF 

We have shown that to accurately describe the hard X-ray 
spectrum of Cep X-4 we require a second absorption feature 
in addition to the prominent CRSF around 30 keV. If this fea¬ 
ture were the fundamental line, it would imply a lower mag¬ 
netic field strength than previously suggested. That is, how¬ 
ever, unlikely as we expect the ratio of the fundamental and 
first harmonic line energy to be close to 2, while we measure 

1.56. Other sources also show deviations from the expected 
factor, however, they are typically much smal ler (< 10%, e.g., 
[Pottschmidt et al.]|2005[ IMiiller et al.||2013] l. The deviation 
we observe is to o large to be explained by relativistic effects 
( |Meszaros|1992| i. 

The second feature is also unlikely to result from sampling 
different accretion columns or different regions. If that were 
the case, the strength of the 30keV line would be difficult to 
explain as this scenario implies that the line should only be 
present at some phase intervals. In order to further investigate 
this we will present results of phase-resolved spectroscopy in 
a separate paper (Bhalerao et al., in prep.). The production 
of two CRSFs of such different energies in different accretion 
columns would also indicate a strong deviation fro m a simple 
dipole magnetic field, which is not expected (see |Nishimura 
2005 1 for a discussion of the influence of multipole helds on 
tHeT^SF). 

The most likely explanation for the second absorption fea¬ 
ture is a deviation of the shape of the CRSF from a smooth 
Gaussian or Lorentzian profile. This deviation is highlighted 
in Figure where we plot the data-to-model ratio using the 
best-fit model with both lines removed. We superimpose 
the ratio between this model and the best-fit model with the 
secondary feature removed, which implicitly shows a simple 
symmetric line profile. From this ratio the data clearly deviate 
at the red wing. A similar deviation might be present on the 
blue side, however, the data quality does not us allow to con¬ 
strain this. Adding a feature there with a similar optical depth 
does not change the statistical quality of the fit significantly. 

Prominent emission wings would be the most obvious 
explanation for deviations from a smooth line. However, 


Fig. 3.— Residuals of observation 1 in terms of data-to-model ratio. The 
ratio was calculated by setting the strength of both absorption features to 
zero. FPMA data are shown in red, FPMB in blue. The best-fit model is 
shown in gray. The green line shows the ratio when only setting the second 
absorption feature to zero, i.e., when assuming a symmetric line profile. The 
data were strongly rebinned for the plot. 

such emission wings typically require a harder spectrum than 
we observe to spawn en ough photons to become significant 
(SchonheiT et al. 2007| l. On the other hand, depending on 
geometry and optical depth, photon spawning can lead to dis¬ 
torted shapes of the scatt ering trough wit hout creating mea¬ 
surable emission wings ( [Schwann 2010| l. As can be seen 
in the inset of Figure and in Figure the second absorp¬ 
tion feature is located exactly at the energy where the primary 
CRSF starts to produce significant deviations from the contin¬ 
uum spectrum, supporting this interpretation. 


4.2. Luminosity dependence of the CRSF 

Thanks to NuSTAR ’s energy resolution, we are able to con¬ 
strain the centroid of the line energy to better than +0.3 keV. 
This allowed us to measure a significant change between the 
two observations, which appears to be correlated with lumi¬ 
nosity. To highlight that correlation and put it into context, 
we show the CRSF energy as a function of luminosity in Fig- 
urefflfor different sources. The energy decreases with declin¬ 
ing mix, leading to a positive correlation between luminosity 
and CRSF energy. The luminosity of Cep X-4 is between 
(l-6)xl0^® erg s“' for a distance of 3.8kpc, which is below 
the theoretical l imit for the formati on of a shock in the ac¬ 
cretion column ( Becker et al.||2012 l. We therefore expect a 
constant CRSF energ y uncorrelated with lu minosit y, as ob¬ 
served in A 0535+26 ^Caballero et al.|2007| but see [Sartor^ 
Jourdain & Roques|2015 for indications of a correlation). In 
recent NuSTAR observations of Vela X-1, [Fiirst et aT] ( [2014a[ l 
show that this source also has a significant positive correlation 
between CRSF energy and luminosity at luminosities below 
10^^ erg s“^. 

[Fiirst et al.[ ( [2014a| l follow the calculations of [Becker et al.[ 
( [2012[ ) but allow for wind-accretion instead of disk-accretion 
and assume a massive neutron star (around 2 Mq in Vela X-1). 
This results in a nar rower accretion column than assumed in 
Becker et al.[ ([2012[l, which in turn decreases the luminosity 
threshold required for shock formation significantly, moving 
the measured luminosities of Vela X-1 partially above it (see 
the dot-dashed line in Figure [^. Cep X-4 also crosses this 
adopted line, although the mass of the neutron star in this sys¬ 
tem is not known. Additionally, Cep X-4 is a Be-system, and 
so we expect that the accreted matter forms a temporary ac¬ 
cretion disk around the neutron star, which leads to a differ- 

































































Distorted cyclotron line profile in Cep X-4 
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Fig. 4.— Dependence of the CRSF energy on luminosity for different 
sources. Cep X-4 data are plotted as red diamonds with a best-fit linear cor- 
relation superimpo sed (dashed line). Vela X-1 data (green asterisk) are from 
|Furst et al.|l2014al and show the energy of t he first harmonic CRSF divided 
by 2. AOS35-t-26 aata (blue crosses) are from lCab allero et S112007 1, H er X- 
1 data (magenta down-pointing triangles) are rrom[tteubertetalJ|20O/Land 
GX 304-1 (orange up pointing triangles) are from| Yamamoto et al.^zOllJ 
and|Klochkov et al.|(2012). The dotted line indicates the original Coulomb 
stopping luminosity presented by |Becke r et al.H2 0121 and the dot-dashed line 
indicates the updated values for Vela X-1 ^Ftirst et al.|2014^. 


ent a ccretion geometry than in purely wind accreting systems 
(e.g. Ghosh & Lamb 1979[ Okazaki, Hayasaki & Moritani] 
|2013 I. If the accretion geometry could be constrained, tins 
correlation might indicate that the neutron star in Cep X-4 is 
also massive, with M > 2Mo. 


4.3. Summary 

NuSTAR has revealed two new interesting features of the 
CRSF in Cep X-4: a distorted profile and a luminosity de¬ 
pendence of the line’s energy. This makes Cep X-4 the first 
system where a significant deviation from a symmetric line 
profile has been measured in the phase averaged spectrum 
and the second system where a positive correlation between 
CRSF energy and luminosity has been found at luminosities 
below 10^^ erg s ’. The latter discovery challenges the cur¬ 
rent understanding of the accretion column geometry, as the 
line forming region is expected to be at the neutron star sur¬ 
face and therefore independent of luminosity. 

The discovery of a complex line profile on the other hand 
is a good qualitative confirmation of theoretically predicted 
line-profiles when taking p hoton-spawning and magnetic field 
gradie nts into account (e.g. Nishimura|2005 |Schonherr et al.| 
2007|l. By combining detailed calculations of the line prohle 
with sophisticated light bending calculations, it should be pos¬ 
sible to improve our understanding of the emission geometry 
in Cep X-4 and other neutron star systems. 
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